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Abstract 

A macroscopic theory for closure relations, referred to as the thermodynamic field theory (TFT), 
applied to systems out of Onsager's region, has been proposed a decade ago. The aim of this 
theory was to determine the nonlinear flux-force relations that respect the thermodynamic the- 
orems for systems far from equilibrium. A new formulation of the TFT is presented herein. In 
this new version, one of the basic restrictions in the old theory, namely a closed-form solution for 
the skew-symmetric piece of the transport coefficients, has been removed. In addition, the general 
covariance principle is replaced by the De Donder-Prigogine thermodynamic covariance principle 
(TCP). The introduction of TCP requires the application of an appropriate mathematical formal- 
ism, which is referred to as the iso-entropic formalism. The validity of the Glansdorff-Prigogine 
Universal Criterion of Evolution, via geometrical arguments, is proven. A new set of closure equa- 
tions determining the nonlinear corrections to the linear ("Onsager") transport coefficients is also 
derived. The geometry of the thermodynamic space is non-Riemannian tending to be Riemannian 
for hight values of the entropy production. In this limit, we obtain again the transport equations 
found by the old theory. Applications of the theory, such as transport in magnetically confined 
plasmas, materials submitted to temperature and electric potential gradients or to unimolecular 
triangular chemical reactions can be found at references cited herein. Transport processes in toka- 
mak plasmas are of particular interest. In this case, even in absence of turbulence, the state of the 
plasma remains close to (but, it is not in) a state of local equilibrium. This prevents the transport 
relations from being linear. 
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I. INTRODUCTION 



It is well known that the basic theory of dynamical systems should provide with an al- 
gorithm for the determination of the moments of the particle distribution functions /" (i.e., 
the average values of the power of particle momenta p), which are determined by the (fluc- 
tuating) fields through the kinetic equations. In the case of turbulent plasmas, for example, 
the most fundamental approach is the study of the stochastic kinetic equation coupled to 
the stochastic Maxwell equations. Such a self-consistent theory should not require any arbi- 
trary assumption: it should produce equations of evolution for all the moments. In practice, 
however, an exact solution of this problem is impossible. Indeed, the equations of evolution 
of the moments have a hierarchical structure: the determination of a moment of order n 
requires the knowledge of order n + 1. Hence, the equations for the third moments will 
involve the fourth moments, and so on ad infinitum. Because of these difficulties, the fun- 
damental studies, in spite of their basic importance, can not easily produce explicit results 
that can be directly compared to experiments. In order to obtain such results, one is led 
to make compromises: we must introduce additional simplifying assumptions allowing to 
truncate the hierarchy. As a result, we obtain a set of dynamical moments equations with 
a number of undetermined quantities: the equations are not closed. These quantities are 
of four kinds: thermodynamic quantities (such as temperature, pressure etc.), electromag- 
netic fields, moments-and energy-exchanges (such as the collisional friction forces or the 
collisional particles heat exchange) and fiuxes (such as, the particle fiux, the heat fiux etc.). 
The dynamics of a thermodynamic system is finally based on the set of balance equations 
coupled to a (macroscopic) theory for the closure relations. Thus, in a macroscopic pic- 
ture of thermodynamic systems, the formulation of a theory for the closure relations plays 
a fundamental role. The connection between the macroscopic equation and a microscopic 
distribution of particles should be established analyzing case by case (for example, for mag- 
netically confined plasmas, see ref. [l] and section [Til - subsection The Nonlinear Closure 
Equations) . 

The most important closure relations are the so-called transport equations, relating the 
dissipative fiuxes to the thermodynamic forces that produce them. The latter are related to 
the spatial inhomogeneity, and are expressed as gradients of the thermodynamic quantities. 
The study of these relations is the object of non-equilibrium thermodynamics. Close to 
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equilibrium, the transport equations of a thermodynamic system are provided by the well- 
known Onsager theory. Indicating with and the thermodynamic forces and fluxes, 
respectively, the Onsager relations read 

Jfi = To^j.uX" (1) 

where tq^,^ are the transport coefficients. We suppose that all quantities involved in Eqs ([1]) 
are written in dimensionless form. Matrix tq^i, can be decomposed into a sum of two ma- 
trices, one symmetric and the other skew- symmetric, which we denote with L^j, and /o/^z^, 
respectively. The second principle of thermodynamics imposes that L^j^^ be a positive def- 
inite matrix. In this equation, as in the remainder of this paper, the Einstein summation 
convention on the repeated indexes is adopted. The most important property of Eqs ([ID is 
that near equilibrium, the coefficients t^^ are independent of the thermodynamic forces, so 
that 

The region where Eqs ([2]) hold, is called Onsager 's region or, the linear region. A well- 
founded microscopic explanation on the validity of the linear phenomenological laws was 
developed by Onsager in 1931 [2]. Onsager's theory is based on three assumptions: 1) The 
probability distribution function for the fluctuations of thermodynamic quantities (Temper- 
ature, pressure, degree of advancement of a chemical reaction etc.) is a Maxwellian 2) 
Fluctuations decay according to a linear law and 3) The principle of the detailed balance 
(or the microscopic reversibility) is satisfied. Onsager showed the equivalence of Eqs ([1]) 
and ([2]) with the assumptions l)-3) [assumption 3) allows deriving the reciprocity relations 
T'Oij.u = 'Tou^J\■ Many important theorems have been demonstrated for thermodynamic sys- 
tems in the linear region. Among them, the most irnportant one is the Minimum Entropy 



Production Theorem., showed by Prigogine in 1947 
the Onsager region, for a — a processes 



3|. This theorem establishes that, in 
a thermodynamic system relaxes towards a 
steady-state in such a way that the rate of the entropy production is negative 

da fda , , \ / \ 

< y— = at the steady statej (3) 



where a = L^^X^X'^ indicates the entropy production and t is time. In 1954, Glansdorff 
and Prigogine demonstrated a more general theorem, valid also when the system is out of 
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Onsager's region They showed that, regardless of the type of processes, a thermodynamic 
system relaxes towards a steady-state in such a way that the following quantity V is negative 

V = -^M-^ < {V = at the steady state j (4) 

Inequality (jl]) reduces to inequahty ([3]) for a — a processes in the Onsager region. For 
spatially-extended systems, the expression in Eqs. (jl]) should be replaced by 

r dX^ / \ 

V = j dv < [V = at the steady statej (5) 

where dv is the (spatial) volume element and the integration is over the entire space Q 
occupied by the system in question. J'^{r,t) and A'^(r,t) denote the space-time dependent 
fluxes and forces, respectively. The phenomenological equations are not needed for deriving 
this more general theorem and no restrictions are imposed to the transport coefficients (apart 
from the validity of the second principle of thermodynamics). Therefore, no use is made of 
the Onsager reciprocal relations, nor it is necessary to assume that the phenomenological 
coefficients are constants. The inequality expressed in (jl]) [or in (j5])] is referred to as the 
Universal Criterion of Evolution and it is the most general result obtained up to now in 
thermodynamics of irreversible processes. Out of Onsager's region, the transport coefficients 
may depend on the thermodynamic forces and Eqs (j2]) may loose their validity. This happens 
when the first end/or the second assumption of the Onsager theory [i.e., the above-mentioned 
assumption 1) end/or assumption 2)] are/is not satisfied. Magnetically confined tokamak 
plasmas are a typical example of thermodynamic systems out of Onsager's region. In this 
case, even in absence of turbulence, the local distribution functions of species (electrons 
and ions) deviate from the (local) Maxwellian. After a short transition time, the plasma 



ocal equilibrium; fiuctuations do not decay 
6(1 and section [III - subsection The Nonlinear 



remains close to (but, it is not in) a state o: 
according to a linear law (see, for example. 
Closure Equations). 

Transport, in the nonlinear region, has been largely studied, both experimentally and 
theoretically. In particular, many theories, based on the Fourier expansion of the transport 
coefficients in terms of the thermodynamic forces, have been proposed (see, for example, 

n n n 

refs [7[, [8(1 and [9|). The theoretical predictions are however in disagreement with the 
experiments and this is mainly due to the fact that, in the series expansion, the terms of 
superior order are greater than those of inferior order. Therefore truncation of the series at 
some order is not mathematically justified. 



A thermodynamic field theory (TFT) has been developed in 1999 for proposing a closure 
relations theory for thermodynamic systems out of the Onsager region. In particular, the 
main objective of this work was to determine how the linear flux-force relations [i.e., Eqs ([T])] 
should be "deformed" in such a way that the thermodynamic theorems for systems far from 
equilibrium are respected [l^. The Onsager coefficients enter in the theory as an input in the 
equations and they have to be calculated by kinetic theory. Attempts to derive a generally 



covariant thermodynamic field theory (GTFT) can be found in refs The characteristic 
feature of the TFT is its purely macroscopic nature. This does not mean a formulation based 
on the macroscopic evolution equations, but rather a purely thermodynamic formulation 
starting solely from the entropy production and from the transport equations. The latter 
provide the possibility of defining an abstract space (the thermodynamic space), covered by 
the n independent thermodynamic forces X^, whose metric is identified with the symmetric 
part of the transport matrix. The law of evolution is not the dynamical law of particle 
motion, or the set of two-fiuid macroscopic equations of plasma dynamics. The evolution 
in the thermodynamic space is rather determined by postulating three purely geometrical 
principles: the Shortest path principle, the Skew-symmetric piece of the transport coefficients 
in closed form and the Principle of least action. From theses principles, a set of closure 
equations, constraints, and boundary conditions are derived. These equations determine 
the nonlinear corrections to the linear ("Onsager") transport coefficients. However, the 
formulation of the thermodynamic field theory, as reported in refs lOj, raises the following 
fundamental objection 

There are no strong experimental evidences supporting the requirement that the skew- 
symmetric piece of the transport coefficients is in a closed form. 

Moreover, the principle of general covariance, which in refs jjy] has been assumed to be 
valid for general transformations in the space of thermodynamic configurations, is, in re- 
ality, respected only by a very limited class of thermodynamic processes. In this paper, 
through an appropriate mathematical formalism, the is o-entropic formalism, the entire TFT 
is re-formulated removing the assumptions regarding the closed-form of the skew-symmetric 
piece of the transport coefficients and the general covariance principle (GCP). The GCP 
is replaced by the thermodynamic covariance principle (TCP), or the De Donder-Prigogine 
statement establishing that thermodynamic systems, obtained by a transformation 
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of forces and fluxes in such a way that the entropy production remains unaltered, are ther- 
modynamically equivalent. This principle applies to transformations in the thermodynamic 
space, and they are referred to as the thermodynamic coordinate transformations (TCT). It 
is worthwhile mentioning that the TCP is actually largely used in a wide variety of thermo- 
dynamic processes ranging from non equilibrium chemical reactions to transport processes 
in tokamak plasmas (see, for examples, the papers and books cited in refs and Q). 
To the author knowledge, the validity of the thermodynamic covariance principle has been 
verified empirically without exception in physics until now. 

The analysis starts from the following observation. Consider a relaxation process of a ther- 
modynamic system in the Onsager region. If the system relaxes towards a steady-state along 
the shortest path in the thermodynamic space, then the Universal Criterion of Evolution is 
automatically satisfied. Indeed, in this case, we can write 

J^X^ = {L,, + MX'^X^ (6) 

where the dot over the variables indicates the derivative with respect to the arc parameter 
defined as 

d^^ = L^^dX^dX" (7) 

Parameter <^ can be chosen in such a way that it vanishes when the system begins to evolve 
and it assumes the value, say /, when the system reaches the steady-state. In the Onsager 
region, the thermodynamic space is an Euclidean space with metric L^^,. The equation of 
the shortest path reads = 0, with solution of the form 

X'^ = a\ + 6^ (8) 

where a'^ and b'^ are arbitrary constant independent of the arc parameter. Inserting Eq. ([8]) 
into Eq. and observing that L^^a'^a'^ = 1 and f^^uO'^a'^ = 0, we find 

J^X^ = ^ + T^^^a^'h" (9) 

At the steady state (i.e. for <, = I) J^X^ \st.state= (because V \st.state= 0). Eq. ([9]) can then 
be written as 

P = -(/ - ^) < (with P = J^X^) (10) 

or 



The equation for the dissipative quantity P, when the thermodynamic system relaxes in the 
hnear region, is thus given by Eq. ([9]): 

Also note that a = 2P < i.e., the minimum entropy production theorem is also satis- 
fied during relaxation. Now, our question is: "How can we "deform" the linear flux-force 
relations in such a way that the Universal Criterion of Evolution remains automatically sat- 
isfied, without imposing any restrictions to the transport coefficients, also out of Onsager's 
region . Outside the linear region, one may be tempted to construct a Riemannian space 
(of 3 or more dimensions) which is projectively fiat i.e., having a vanishing Weyl's pro- 
jective curvature tensor. In this case, indeed, there exists a coordinate system such that 
the equations of the shortest path are linear in the coordinates [i.e., the shortest paths are 
given by equations of the form ([8])]. In this respect, we have the following Weyl theorem 



151]: a necessary and sufficient condition that a Riemannian space be projectively fiat is 
that its Riemannian curvature be constant everywhere. On the other hand, to re-obtain 
the Onsager relations, we should also require that, near equilibrium, the Riemannian space 
reduces to a fiat space (which has zero Riemannian curvature). The Weyl theorem can be 
conciliated with our request only if there exists a coordinate system such that Eqs ([2]) are 
valid everywhere, which is in contrast with experiments. Thus one wants the Universal Cri- 
terion of Evolution satisfied also out of the Onsager region, without imposing a priori any 
restrictions on transport coefficients, a non- Riemannian thermodynamic space is required. 
Clearly, a transport theory without a knowledge of microscopic dynamical laws can not be 
developed. Transport theory is only but an aspect of non-equilibrium statistical mechanics, 
which provides the link between micro and macro-levels. This link appears indirectly in 
the "unperturbed" matrices, i.e. the L^^ and the fo^u coefficients, used as an input in the 
equations. These coefficients, which depend on the specific material under consideration, 
have to be calculated in the usual way by kinetic theory. 

In section |Tll we introduce a non- Riemannian space whose geometry is constructed in such 
a way that 

A. The theorems valid when a generic thermodynamic system relaxes out of equilibrium 
are satisfied; 
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B. The nonlinear closure equations are covariant under the thermodynamic coordinate 
transformations (TCT). 

We shall see that the properties of geometry do not depend on the shortest paths but upon a 
particular expression of the affine connection. Our geometry is then of affine type and not of 
projective type. At the end of section [TTl we derive the nonlinear closure equations through 
an appropriate mathematical formalism: the iso-entropic formalism. This formalism allows 
to respect the De Donder-Prigogine statement. New geometrical objects like thermodynamic 
covariant differentiation or the thermodynamic curvature are also introduced. We shall see 
that under the weak-field approximation and when o" ^ 1, but only in these limits, the new 
nonlinear closure equations reduce to the ones obtained in refs llyl. So that, all results found 

n n 

in refs for magnetically confined plasmas, and in refs [16j, for the nonlinear thermoelectric 



effect and the unimolecular triangular reaction, remain valid. In section UTT] it is shown that 
this formalism is able to verify the thermodynamic theorems (in particular, the Universal 
Criterion of Evolution) for systems relaxing out of the Onsager region. Mathematical details 
and demonstrations of the theorems are reported in the annexes. 

II. THE ISO-ENTROPIC FORMALISM 

Consider a thermodynamic system driven out from equilibrium by a set of n independent 
thermodynamic forces {X'^} [fi = 1, ■ ■ ■ n). It is also assumed that the system is submitted 
to time- independent boundary conditions. The set of conjugate flows, {Jfi}, is coupled to 
the thermodynamic forces through the relation 

Jf. = r^uX-" (13) 

where r^,y denote the transport coefficients. The symmetric piece of r^u is denoted with g^^, 
and the skew-symmetric piece as f^^,: 



where 



9f^u = ^{r^u + r^f,) = Quf, (15) 



It is assumed that (yf^i, is a positive definite matrix. With the elements of the transport coef- 
ficients two objects are constructed: operators, which may act on thermodynamic tensorial 
objects and thermodynamic tensorial objects, which under coordinate (forces) transforma- 
tions, obey to well specified transformation rules. 

Operators 

Two operators are introduced, the entropy production operator cr(X) and the dissipative 
quantity operator P{X), acting on the thermodynamic forces in the following manner 

a{X) a{X) = XgX^ 

P{X) :^ P{X) = XrX^ (17) 

In Eqs ( !T7|) . the transport coefficients are then considered as elements of the two n x n 
matrices, r and g. The positive definiteness of the matrix g^^i, ensures the validity of the 
second principle of thermodynamics: a > 0. These matrices multiply the thermodynamic 
forces X expressed as n x 1 column matrices. The dot symbol stands for derivative with 
respect to parameter q, defined in Eq. (1291) . We anticipate that this parameter is invariant 
under the thermodynamic coordinate transformations. Thermodynamic states Xg such that 

P{Xg) = (18) 

are referred to as steady- states. Of course, the steady-states should be invariant expressions 
under the thermodynamic coordinate transformations. Eqs f|T7j) should not be interpreted as 
the metric tensor g^^,, which acts on the coordinates. The metric tensor acts only on elements 
of the tangent space (like dX'^, see the forthcoming paragraphs) or on the thermodynamic 
tensorial objects. 

Transformation Rules of Entropy Production, Forces and Flows 

According to the De Donder-Prigogine statement |l2|, |l^ thermodynamic systems are 
thermodynamically equivalent if, under transformation of fluxes and forces the bilinear form 



of the entropy production, a, remains unaltered 29| . In mathematical terms, this implies: 



a = J^X'' = j;X'^ (19) 
This condition requires that the transformed thermodynamic forces and flows satisfy the 
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relation 

■r. - ax^^^ (20) 

These transformations are referred to as Thermodynamic Coordinate Transformations 
(TCT). The expression of entropy production becomes accordingly 

a = J^X^ = T,,X^^X^ = g^,X>^X- = g'^^X'^X'^ = a' (21) 



9xk - ^^"''^xo:'gx^ ^22) 



From Eqs (l20l) and (1211) we derive 



Moreover, inserting Eqs fl20|) and Eq. fl22l) into relation = {g^^^ + f^^)X'^, we obtain 

r-(a' +f ^^'^^^W- (23) 
J\ - [9xk + Ji^u )X [26) 

or 

f)X^ BX'^ 

A = {9L + fL)X'^ with /L = /..^^ (24) 



Hence, the transport coefficients transform like a thermodynamic tensor of second order 
Properties of the TCT 

By direct inspection, it is easy to verify that the general solutions of equations fl2U]) are 

-2 v3 J^n 



3Q| 



(25) 



where is an arbitrary function of variables X^ /X^~'^ with (j = 2, . . . , n). Hence, the TCT 
may be highly nonlinear coordinate transformations but, in the Onsager region, we may (or 
we must) require that they have to reduce to 

X'^" = c'^X'' (26) 

where are constant coefficients (i.e., independent of the thermodynamic forces). Note 
that from Eq. (l20i) . the following important identities are derived 

d'^X'^ , d'^X^ 

= ° ' ^^^^ 
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Moreover 

dX'^ = ^^dX" 
dX'^ 

dX'i' dX't'dX'^ ^ ' 

i.e., dX^ and d/dX^ transform like a thermodynamic contra-variant and a thermodynamic 
covariant vector, respectively. According to Eq. (1281) . thermodynamic vectors dX^^ define the 
tangent space to Ts. It also follows that the operator P{X), i.e. the dissipation quantity, 
and in particular the definition of steady-states, are invariant under TCT. Parameter 
defined as 

dq^ = g^^dX^dX" (29) 
is a scalar under TCT. The operator O 

O = Xf"-— = X'^— — = O' (30) 

is also invariant under TCT. This operator plays an important role in the formalism. 

Thermodynamic Space, Thermodynamic Covariant Derivatives and Thermody- 
namic Curvature 

A non-Riemannian space with a linear connection F^^ is now introduced. Consider an 
n-space in which the set of quantities is assigned as functions of the n independent 
thermodynamic forces X^, chosen as coordinate system. Under a coordinate (forces) trans- 
formation, it is required that the functions F^^ transform according to the law 



_ dX'^^ dX' dx- dX'^^ d^x-^ 

a/3 A« gj^^ Q^,^ + dX"^dX'P ^ ' 



With the linear connection F^^, the absolute derivative of a thermodynamic contra- variant 
vector along a curve can be defined as 

It is easily checked that, if the parameter along the curve is changed from q to then 
the absolute derivative of a thermodynamic tensor field with respect to q is dq/dg times 
the absolute derivative with respect to q. The absolute derivative of any contra-variant 
thermodynamic tensor may be easily obtained generalizing Eq. (l32l) . In addition, the linear 
connection is submitted to the following basic postulates: 
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1. The absolute derivative of a thermodynamic contra-variant tensor is a thermodynamic 

tensor of the same order and type. 

2. The absolute derivative of an outer product of thermodynamic tensors, is given, in terms 

of factors, by the usual rule for differentiating a product. 

3. The absolute derivative of the sum of thermodynamic tensors of the same type is equal 

to the sum of the absolute derivatives of the thermodynamic tensors. 

In a space with a linear connection, we can introduce the notion of the shortest path defined 
as a curve such that a thermodynamic vector, initially tangent to the curve and propagated 
parallelly along it, remains tangent to the curve at all points. By a suitable choice of the 
parameter the differential equation for the shortest path is simplified reducing to 

+ dg dg "° ^^^^ 

To respect the general requirement A. (see section it is required that the absolute 
derivative of the entropy production satisfies the equality 

More in general, it is required that the operations of contraction and absolute differentiation 
commute for all thermodynamic vectors. As a consequence, the considered space should be 
a space with a single connection. The absolute derivative of a covariant thermodynamic 
vector T^ is then defined as 

ST,_dT, dXP 

The absolute derivative of the most general contra- variant, covariant and mixed thermody- 
namic tensors may be obtained generalizing Eqs (!32l) and (1351) . The derivatives, covariant 
under TCT, of thermodynamic vectors, are defined as 

dTf" 
dT 

= ^ - r^.T„ (36) 
For the entropy production, it is also required that 

CTlMli. = Cr|i.|M (37) 
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More in general, Eq. (137|) should be verified for any thermodynamic scalar T. This postulate 
requires that the linear single connection F^^ is also symmetric i.e., F^^ = F|^^. A non- 
Riemannian geometry can now be constructed out of n^(?2 + l)/2 quantities, the components 
of F^^, according to the general requirements A and B mentioned in the introduction. 
In the forthcoming paragraph, the expression of the affine connection F^^ is determined from 
assumption A. In section IIIII it is shown that the Universal Criterion of Evolution, applied 
to thermodynamic systems relaxing towards a steady-state, is automatically satisfied along 
the shortest path if, in case of symmetric processes (i.e., for a — a processes), we impose 

[dX^ + dX^-dX^)+^^^^9a,) (38) 
where 0{g^,) ^ 



In the general case, we have 



ATM' 

4 



+ (39) 



2a dX^^J 
where the thermodynamic Christoffel symbols of the second kind are introduced 

2^ \dX^ dX^ dx^i ^ ' 



a(3 



and matrix N^'^ is defined as 



N^,=g^, + -f^^X^X, + -U^X^X^ with iV^'^: iV'^^AT,, = (41) 
a a 

In appendix |A] it is proven that the affine connections Eqs (138!) and (139|) transform, under a 
TCT, as in Eq. ( 13T|) and satisfy the postulates 1., 2. and 3. From Eq. ( l4Ti) we easily check 
that 

N^u = N,^ 

NfiuX" = {^^y + —ffj^f^X'^Xi, + —fy^X'^Xf^X" = Jfj. 

N^,X^ = K^X^ = J, (42) 
N^.X'X^' = J^X^ = a 
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While 



Nf^'^J^ = Nf^'^N^^X" = N'^'^N^f^X" = X'^ (43) 
N^^JJ^ = X^J^ = a 

The shortest path is the same for two symmetric connections whose coefficients are related 

as 

= + ^'a^P + ^p^a (44) 

where ipa is an arbitrary covariant thermodynamic vector and 5^ denotes the Kronecker 
tensor. In literature, modifications of the connection similar to Eqs (1441) are referred to as 
projective transformations of the connection and tpa the projective covariant vector. The 
introduction of the affine connection gives rise to the following difficulty: the Universal 
Criterion of Evolution is satisfied for every shortest path constructed with affine connections 
r^^, linked to F^^ by projective transformations. This leads to an indetermination of the 
expression for the affine connection, which is not possible to remove by using the De Donder- 
Prigogine statement and the thermodynamic theorems alone. This problem can be solved by 
postulating that the nonlinear closure equations (i.e., the equations for the affine connection 
and the transport coefficients) be symmetric and projective-invariant (i.e., invariant under 
projective transformations). 

For any covariant thermodynamic vector field T^, we can form the thermodynamic tensor 



R'^Xk the following manner 17 1 



T.|A|. - n\.\x = T,R^^,^ (45) 

where 



dX^ dX'^ ~'~ ^ ^''^ ^ '^'^'^ ^'^ 
with R'^xk satisfying the following identities 



^uXk ~ ^ukX 



Kx. + RLu + Ku. = (47) 



I +R^ n+R'' = 

uXkIt] UKri\X uriX\K 
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By contraction, we obtain two distinct thermodynamic tensors of second order 

-^!^A J^^Xf, Q-^^ ^ J- u^l^ nx J- i^x^ w 

with Fxu being skew-symmetric and R^x asymmetric. Tensor R^x can be re-written as 
R.x = Byx + Fxu where 

r> _ r> _ _ J >Hi \ _ uX , pr? _ p*? p/^ rAQ\ 

J^vx - oxu - 2\ Qx^ dX" J dXt" ''^ ^ ^ 



Hence, Fxu is the skew-symmetric part of R^x |3l|- It is argued that the closure equations 
can be derived by variation of a stationary action, which involves R^x- Symmetric and 
projective-invariant closure equations may be obtained by adopting the following strategy: 
1) a suitable projective transformation of the affine connection is derived so that R^x be 
symmetric and Fxu be a zero thermodynamic tensor and 2) the most general projective 
transformation that leave unaltered Riyx and Fxu (= 0) is determined. By a projective 
transformation, it is found that 

B,x = B,x + ^(Ifl - ^-^a) - - ^.^a) 



^- = ^- + — laf^-af^J (50) 



Eq. fHS]) shows that Fxv can be written as the curl of the vector aj,/2 defined as |20| 



= TL - <( > (51) 



Consequently, by choosing 




^^ = -;^(rL-<! )] (52) 



we have Fxu = and R^x = B^x- From Eqs fl50l) . we also have that the thermodynamic 
tensor Ri^x remains symmetric for projective transformations of connection if, and only if, 
the projective covariant vector is the gradient of an arbitrary function of the X's |2o|]- In 
this case, the thermodynamic tensor Fxv remains unaltered i.e., Fxu = 0. Hence, at this 
stage, the expression of the affine connection is determined up to the gradient of a function, 
say 0, of the thermodynamic forces, which is also scalar under TCT. Let us impose now the 
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projective-invariance. Eqs (150|) indicate that a necessary and sufficient condition that R,y\ 

be projective-invariant is that 

d'^(j) d(J) d(j) ^ 

U with 



dX^dX'^ dX^ dX'^ 


= (in the Onsager region) 



(53) 







where is a function, invariant under TCT. The solution of Eq. (jnSD is = everywhere. 
The final expression of the affine connection for symmetric processes reads then 



a/3 1 2a 2{n + l)cr 



5^X''0{gp,)+5>;X''0{g^,) 



(54) 



The general case is given by 



A I N^''^ 
«A s > + ^X^Oig^fs) + 



2a \dXP^ dX-) 



(55) 



where 



n+1 



A 



2(n + l)a^^'^^^ 2(n + l)a " \dX^ ^ dX^ 



-f..x^xH^ + ' 



(56) 



2(n + l)a''"' KdX'^ dXn ) n + l dX'^ 

Note that the thermodynamic space tends to reduce to a (thermodynamic) Riemannian 
space when a^^ ^ 1. The following definitions are adopted: 

• The space, covered by n independent thermodynamic forces X'*, with metric tensor g^y 
and a linear single connection given by Eq. fl55l) . may be referred to as thermodynamic 
space Ts (or, space of the thermodynamic forces). 



In Ts, the length of an arc is defined by the formula 

L = 



gf,u— — ) 



The positive definiteness of matrix (yf^j, ensures that L > 0. Consider a coordinate system 
X^^, defining the thermodynamic space Ts. 
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• All thermodynamic spaces obtained from Ts by a TCT transformation, may be called 

iso-entropic spaces. 

In the TFT description, a thermodynamic configuration corresponds to a point in the ther- 
modynamic space Ts. The equilibrium state is the origin of the axes. Consider a thermody- 
namic system out of equilibrium, represented by a certain point, say a, in the thermodynamic 
space 

• A thermodynamic system is said to relax towards another point of the thermodynamic 
space, say b, if it moves from point a to point b following the shortest path (133|) . with 
the affine connection given in Eq. fl55l) . 

• With Eq. (!55|) . Eqs (!36|) may be called the thermodynamic covariant differentiation 
of a thermodynamic vector while Eqs. fl32l) and fl35|) the thermodynamic covariant 
differentiation along a curve of a thermodynamic vector. 

• With affine connection Eq. fl55l) . R'^xk called the thermodynamic curvature 
tensor. 

• The scalar R obtained by contracting the thermodynamic tensor R^x with g"^^ (i.e. 
R = Ruxg"^) may be called the thermodynamic curvature scalar. 

The Principle of Least Action 

From expression (l55l) . the following mixed thermodynamic tensor of third order can be 
constructed 

^a, = N g., <> + -^X^Oig^,) + -^X^X + ^ 

+ + M (58) 



= = (59) 



2a ■"^^ \dXP dX»J 
This thermodynamic tensor satisfies the important identities 

a/3 

Again, from the mixed thermodynamic tensor of fifth order can be constructed 
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By contraction, a thermodynamic tensor of third order, a thermodynamic vector and a 
thermodynamic scalar can be formed as follows 



= 5f = ^Kpg"-' (61) 



S ^ S^^'^l = 2M/L*L^7'^^ 



The following postulate is now introduced: 

There exists a thermodynamic action I, scalar under TCT , which is stationary with respect 
to arbitrary variations in the transport coefficients and the affne connection. 
This action, scalar under TCT, is constructed from the transport coefficients, the affine 
connection and their first derivatives. In addition, it should have linear second derivatives 
of the transport coefficients and it should not contain second (or higher) derivatives of the 
affine connection. We also require that the action is stationary when the affine connection 
takes the expression given in Eq. ( |55l) . The only action satisfying these requirements is 



I = I [R,u - (r^^ - Ti,)Sf^^] g^^^gd-X (62) 

where d^X denotes an infinitesimal volume element in Ts and F^^ is the expression given in 
Eq. ( l55i) i.e., f^^ = \E'^^ + l/Ti^j- To avoid misunderstanding, while it is correct to mention 
that this postulate affirms the possibility of deriving the nonlinear closure equations by a 
variational principle it does not state that the expressions and theorems obtained from the 
solutions of these equations can also be derived by a variational principle. In particular the 
well-known Universal Criterion of Evolution established by Glansdorff-Prigogine can not be 
derived by a variational principle (see also section IIIB . 

The Nonlinear Closure Equations 

The transport coefficients and the affine connection should be considered as independent 
dynamical variables (as opposed to X^, which is a mere variable of integration) 21|. There- 
fore, the action (|62l) is stationary with respect to arbitrary variations in g^j_y, f^^, and F^^,. As 
a first step, we suppose that the transport coefficients and the affine connection be subject to 
infinitesimal variations i.e., g^.^ g/iu + Sgf^u, fi^u f^iu + 5f^u and FJ;^ FJ^^ + 5FJJ^, where 
^Q^iui Sffiu and (5F^^, are arbitrary, except that they are required to vanish as | X^ |— > 00. 
Upon application of the principle of stationary action, the following nonlinear closure equa- 

19 



tions (i.e., the equations for the transport coefficients and the affine connection) are derived 
(see appendix [B|) : 

1 iSV^ 
^9tiun- ^^^^ 

Sf-^ = (63) 



where the variations of the affine connection (]55l) with respect to the transport coefficients 
appear in the ffist two equations. Notice that R^i, — \9^uR does not coincide with Einstein's 
tensor (see also Appendix [C]) . From the ffist equation of Eqs fl63l) . the expression for the 
thermodynamic curvature scalar is given as 

R = -^g^^^sf-^ (64) 
n-2 ^ 6gi"' ^ ' 

The third equation of Eqs (1631) can be re-written as 

9i.v,\ - ^"x9au - Kx9a^, = -'^°x9au " '^uX9a^, (65) 

where the comma (, ) denotes partial differentiation. Adding to this equation the same equa- 
tion with /i and A interchanged and subtracting the same equation with u and A interchanged 
gives 

9^lu,x + 9XU,,, - 9f.x,u = '^9auTx^, - 2^^!.^^ (66) 

or 

n, = {x,} + n, = n, (67) 

Hence, action Eq. ( I62l) is stationary when the affine connection takes the expression given 
in Eq. (155|) . For a — a processes, close to the Onsager region, it holds that 

9^^ly = + h^^ + O(e^) 

A. = 0(e) with .^Maxi'^-^^^^^^^^^ (68) 

where Aa- = 1/cr and h^j^ are small variations with respect to Onsager's coefficients. In this 
region, Eq. ( !62l) is stationary for arbitrary variations of h^y and VJ^^. It can be shown that 



m 



tXk ^ h^jy , tXk 9 hxn _ jXk ^ hxy _ xk 9 hxfj, ^ Q , ^/ 2^ 
dX^dX'^ dXi'dX'' dX'^dX'' dX^dX'' ^ 

_ 1 
- 2" 



9^.u;x = + 0ie') or F)!, = -L''^(/i^,,, + /i,,,^ - V^) + ©(e^) (69) 
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where the semicolon indicates the covariant derivative performed with the Christoffel sym- 
bols. Eqs fl68|) should be solved with the appropriate gauge-choice and boundary conditions. 
The validity of Eqs fl69p has been largely tested by analyzing several symmetric processes, 



such as the thermoelectric effect and the unimolecular triangular chemical reactions [10 1. 
More recently, these equations have been also used to study transport processes in magnet- 
ically confined plasmas. In all examined examples, the theoretical results of the TFT are in 
line with experiments. It is worthwhile mentioning that, for transport processes in tokamak 
plasmas, the predictions of the TFT for radial energy and matter fluxes are much closer 
;o the experimental data than the neoclassical theory, which fails with a factor 10'^ -r- 10^ 

in 

1|, [6[. The physical origin of this failure can be easily understood. As mentioned in the 
introduction, even in absence of turbulence, the state of the plasma is close to, but not in, 
a state of local equilibrium. Indeed, starting from an arbitrary initial state, the collisions 
would tend, if they were alone, to bring the system very quickly to a local equilibrium state. 
But slow processes, i.e. free-flow and electromagnetic processes, prevent the plasma from 
reaching this state. Hence, the probability distribution for the fluctuations of the thermody- 
namic quantities deviates from the Maxwellian preventing the thermodynamic fluxes from 
being linearly connected with the conjugate forces. In tokamak plasmas, the thermodynamic 
forces and the conjugate flows are the generalized frictions and the Hermitian moments, re- 
spectively 1^. In the neoclassical theory, the flux-force relations have been truncated at the 
linear order (ref., for example, to [26|), which is inconsistent with the statement that the 
distribution functions of particles (electrons and ions) are not (locally) Maxwellian. This 
may be one of the main causes of the strong disagreement between the neoclassical previsions 
and the experimental profiles [l|. It is, however, important to mention that it is well 
accepted that another main reason of this discrepancy is attributed to turbulent phenomena 
existing in tokamak plasmas. Fluctuations in plasmas can become unstable and therefore 
amplified, with their nonlinear interaction successively leading the plasma to a state, which 
is far away from equilibrium. In this condition, the transport properties are supposed to 
change significantly and to exhibit qualitative features and properties that could not be 
explained by collisional transport processes, e.g. size-scaling with machine dimensions and 
non-local behaviors that clearly point at turbulence spreading etc. The scope of the work 
cited in ref. [l| is mainly to demonstrate that collisional transport processes in fusion plas- 
mas can be computed via a nonlinear theory on a more rigorous and sound basis than that 
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provided by the well known classical and/or neoclassical theory. The proposed approach 
includes prior known results as a limiting case where nonlinear and non-local effects in coUi- 
sional transport processes can be ignored. More generally, the TFT estimates of collisional 
transport fluxes can be amplified by up to two or three orders of magnitude with respect 
to the classical/neoclassical levels in the electron transport channel, while ions corrections 
are much smaller. However, TFT collisional transport levels remain a fraction of the values 
observed experimentally, confirming that turbulent transport is the generally dominant pro- 
cess determining particle and heat fluxes in magnetically confined plasmas. In this specific 
example, the nonlinear corrections provide with an evaluation of the (parallel) Hermitian 
moments of the electron and ion distribution functions [l| . 

Some Remarks on Spatially-Extended Thermodynamic Systems 

The macroscopic description of thermodynamic systems gives rise to state variables that 
depend continuously on space coordinates. In this case, the thermodynamic forces possess 
an infinity associated to each point of the space coordinates. The system may be subdivided 
into cells {NxNxN in three dimensions), each of which labeled by a wave- number k, and 
we follow their relaxation. Without loss of generality, we consider a thermodynamic system 
confined in a rectangular box with sizes Ix, ly and Iz- Denoting with k the wave-number 



Ht^t^t) ^^^^ 

''y '•z 



^ = 0,±l,---±Ar^ 

= 0,±l,---±iV, (70) 



= 0, ±1, ■ ■ ■ ± Nz 
the fiuxes and forces, developed in (spatial) Fourier's series, read 



N 



^/^(i"' ^) ^Mk) (^) exp(zk ■ r) 

n=-N 
N 

X^{v,t)=Y, Xf,,)(t)exp(zk'-r) (71) 

n'=-N 

where, for brevity, n and N stand for n = {nx,ny,nz) and N = {N^, Ny, N^), respectively. 
The Fourier coefficients are given by 

■Jfi(k)it) = ^ j J^,{v,t) exp{~ik- r)dv 

X^^,)it) = ^J^X^{r,t)expi-tk'-v)dv (72) 
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In particular, the contributions at the thermodynamic hmit (i.e., for k — 0) are expressed 

as 

Jmo)W = ^ / Mr,t)dv = J^{t) 

%){t) = ^ I^X'^{r,t)dv = X^{t) (73) 
The entropy production and the fluxes-forces relations take, respectively, the form 

a(r,t) = J'^(r,t)A'^(r,t) >0 

MT,t) = T^,{r,t)X''{r,t) (74) 

Considering that 

rlx ply rlz 

/ / / exp[i(k + k') ■ r]dv = 5k+k',o with (75) 
Jo Jo Jo 

I if k + k' 7^ 

i^k+k'.o = \ and Q = IxLh 

[l if k + k' = 

from the first equation of Eq. flTll) we also find 

f J,{r,t)A:^ir,t) di; = r](j^(o)(t)X5)(t) + 5^J^(k)(t)Xf'_^)(t)) >0 (76) 

On the other hand, we have 

Jm (t) = V(o) (t)Xfo) (t) + J2 V(k) (i)^(-k) (t) (77) 

where 

V(k)W = ^ / r^u{r,t)exp{-ik- T)dv (78) 
" Jn 

Eq. (176|l can then be brought into the form 

a dv = %,^(o)(t)Xj)(t)X(;))(t) 

+l^$^(v(k)(t)X('^-k)(t)^(o)(i) + ^Mk)(i)^r-k)W) > (79) 



where 



9f,u(k)it) = ^ [ g^^{r,t)exp{-ik-r)dv with 
" Jn 

g^.{r,t) = l[r^,(r,t) + r,^(r,t)] (80) 
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In Eq. (179ll . the first term is the contribution at the thermodynamic limit whereas the second 
expression reflects the interactions between the k-cell and the other cells. In a relaxation 
process, contributions from different wave-numbers are negligible with respect to those with 
same wave-numbers {the slaving principle [22]) and, hence, we flnally obtain 

f advc:^ fi^^,(o) (t)X('^o) (^) > V XJ) (t) (and a ^ 0) (81) 

Last inequality is satisfled for any X^g^ (t) if, and only if 

9^u{o)it) = ^ y Gf,uir,t) dv = g^^{t) (82) 

is a positive definite matrix. Therefore, for spatially-extended thermodynamic systems, we 
replace X^(t) — ^ X^'^j(t) and T^y{t) Tf^u{k)it)- Under these conditions, Eqs fl63|) determine 
the nonlinear corrections to the Onsager coefficients while Eqs (|32l) and Eqs (!36l) . with affine 
connection Eq. ( l55i) . are the thermodynamic covariant differentiation along a curve and the 
thermodynamic covariant differentiation of a thermodynamic vector, respectively. 

The Privileged Thermodynamic Coordinate System 

By deflnition, a thermodynamic coordinate system is a set of coordinates deflned so that 
the expression of the entropy production takes the form of Eq. ( |T9i) . Once a particular set 
of thermodynamic coordinates is determined, the other sets of coordinates are linked to the 
flrst one through a TCT [see Eqs (120|) ]. The simplest way to determine a particular set of 
coordinates is to quote the entropy balance equation 

^ + V-J. = a (83) 

where ps is the local total entropy per unit volume and is the entropy flux. Let us consider, 
as an example, a thermodynamic system conflned in a rectangular box where chemical 
reactions, diffusion of matter, macroscopic motion of the volume element (convection) and 
leat current take place simultaneously. The entropy flux and the entropy production read 



23|, M 



i i 



(85) 
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where /ij, piSi and Ai are the chemical potential, the local entropy and the affinity of species 
"z", respectively. Moreover, Fj indicates the external force per unit mass acting on "i", Iljj 
the components of the dissipative part of the pressure tensor M.^^ {J^ij = p6ij + Uij] p is 
the hydrostatic pressure) and Vj is the component of the hydrodynamic velocity (see, for 
example, ref. 25|). The set of thermodynamic coordinates is given as 



T' T 



(86) 



For this particular example, this set may be referred to as the privileged thermodynamic coor- 
dinate system. Other examples of privileged thermodynamic coordinate system, concerning 



magnetically confined plasmas, can be found in refs 



and 



26|. 



III. THERMODYNAMIC THEOREMS FOR SYSTEMS OUT OF EQUILIBRIUM 



In 1947, Prigogine proved the minimum entropy production theorem [3|, which concerns 
the relaxation of thermodynamic systems near equilibrium. This theorem states that: 

Minimum Entropy Production Theorem (MEPT) 

For a — a processes, a thermodynamic system, near equilibrium, relaxes to a steady-state 
Xs in such a way that the inequality 

'i<-" («^) 

is satisfied throughout the evolution and is only saturated at Xg. 

The minimum entropy production theorem is generally not satisfied far from equilibrium. 
However, P. Glansdorff and I. Prigogine demonstrated in 1954 that a similar theorem con- 
tinues to hold for any relaxation to a steady-state, which reads 5j 

Universal Criterion of Evolution (UCE) 

When the thermodynamic forces and conjugate flows are related by a generic asymmetric 
tensor, regardless of the type of processes, for time-independent boundary conditions a ther- 
modynamic system, even in strong non- equilibrium conditions, relaxes towards a steady-state 
in such a way that the following universal criterion of evolution is satisfied: 

dX^' , , 

P . J,— < (88) 

This inequality is only saturated at Xg. 
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For a — a processes, the UCE reduces to the MEPT in the Onsager region. As mentioned 
in the introduction of this manuscript, Glansdorff and Prigogine demonstrated this the- 
orem using a purely thermodynamical approach. In this section we shall see that if the 
system relaxes towards a steady-state along the shortest path then the Universal Criterion 
of Evolution is automatically satisfied. 

By definition, a necessary and sufficient condition for a curve to be the shortest path is that 

it satisfies the differential equation 

cPX^' „ rfX" dXf^ , , dX^ , , 

+ r::,-3--3- = ^W-3^ (89) 



dt2 dt dt dt 

where ip{t) is a determined function of time. If we define a parameter g by 



= cexp f ip*dt with ^p* = if — 2ipi, ^ (90) 



where c is an arbitrary constant and ipu the projective covariant vector, Eq. ( 1891) reduces 
to Eq. (133!) with P^^ given by Eq. (1391) . Parameter g is not the affine parameter s of the 
shortest path. The relation between these two parameters is 

= bj exp(^-2 J 'ip^dX'^^ds (91) 

where b is an arbitrary constant. Eq. ( l90l) allows us to choose the parameter g in such a 
way that it increases monotonically as the thermodynamic system evolves in time. In this 
case, c is a positive constant and, without loss of generality, we can set c = 1. Parameter g 
can also be chosen so that it vanishes when the thermodynamic system begins to evolve and 
it takes the (positive) value, say /, when the system reaches the steady-state. Multiplying 
Eq. fl5^ with the flows and contracting, we obtain 

'd^^ ^ dg dg 



d'^X" dp fdq\-i dX" dX^ ^^dga,x dX^^dX^ .df^x 



However 

^ dg^ dg \dg) dg dg dX^ dg dg OX^ ^ ^ 

where P = Jfi^^ and after taking into account the identities fiiv^^^^ = and 
Qfiu^^^^ = 1. In addition, recalling Eq. and the relations X^X'^ = a and 

f^yX^X^ = 0, it can be shown that 

dX^dX^ ^ dX-dX^ .dg^ dJC^dX^ .dU .q,^ 
^ dg dg dg dg dXP dg dg dXP ^ ' 
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Summing Eq. (p3|) with Eq. and considering Eq. (p2|) . gives 

dP _ /dqy 
dg \dgJ 

Integrating Eq. (!95l) ] from the initial condition to the steady-state, we find 



(95) 



P{X,)-P= I l^VdgyO (96) 



-dg 

From Eq. (fTSl) we have -P(Xs) = P{Xs)dq/dg = 0, so we finally obtain 

P^.,^^- f(*)\„<0 (97) 



where the inequality is only saturated at the steady-state. Recalling Eq. (pOl) . the inequality 
estabhshed by the UCE can be derived 

Eq. fl95|l can be re-written as 



dgJ 



P 



dq 
dg 



(99) 



This equation generalizes Eq. (fT2|) . which was valid only in the near-equilibrium region 
(notice that, in the linear region, d^jdg = 1/h = const.). Integrating Eq. fl99|) . the expression 
of the dissipative quantity P is derived 

P = -('^\ t(^)d^' = -(g ^^^]'^'[\g ^^^fd^' < (100) 
^d(;) \dg) ^ dg dg ) dg dg ) ^ ~ 

On the right, it is understood that the X's are expressed in terms of g{(^). Eq. f llOOp 

generalizes Eq. ( fTOl) . which was valid only in the linear region. For a — a processes in the 

Onsager region, Eq. fl98|) implies the validity of the inequality (|87|) . Indeed, Eq. flMl) gives 

'jL^^JL^j/^^X^'-^ = 2j/-^^X^X^'-^ (101) 

dg dg dg dg dg dg 

In the linear region, the coefficients of the affine connection vanish. Eq. fllOip is simplified 
reducing to 

da da dg / dX^ dg\ ^ 

2 J^-— -f =2P<0 (102) 



dt dg dt V ^ dg dt . 
where the inequality is saturated only at the steady state. 

Let us now consider spatially extended thermodynamic systems. In space-dependent sys- 
tems, the dissipative quantity should be expressed in the integral form 

V= I J^{r,t)dtX^'{T,t) dv (103) 
Jq 
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where dtX^ = dX^/dt. In terms of wave- vectors k, Eq. fllOSp can easily be brought into the 
form 

V = fi(j^(o)(t)rf*x;;)(t) + ^ j^(k)(t)rfiX(%(t)) (104) 

where Eq. (175|) has been taken into account. As aheady mentioned in section Hit in a 
relaxation process, contributions from different wave-numbers are negligible with respect to 
those with same wave-numbers [22]. Thus, considering Eqs (|75]) . we finally obtain 

P = [ J^iv, t)dtX^'{v, t) dv - nj^{t)dtX^{t) < (105) 
Jn 

where inequality (198|1 has been taken into account. It is therefore proven that the Universal 
Criterion of Evolution is automatically satisfied if the system relaxes along the shortest path. 
Indeed it would be more exact to say: the affine connection, given in Eq. (!39|) . has been 
constructed in such a way that the UCE is satisfied without imposing any restrictions on the 
transport coefficients (i.e., on matrices g^^, and f^^). In addition, analogously to Christoffel's 
symbols, the elements of the new affine connection have been constructed from matrices g^^^ 
and /^jy and their first derivatives in such a way that all coefficients vanish in the Onsager 
region. Eq. (139!) provides the simplest expression satisfying these requirements. 

The Minimum Rate of Dissipation Principle (MRDP) 
In ref . jll| the validity of the following theorem is shown: 
The generally covariant part of the Glansdorff-Prigogine quantity is always negative and 
is locally minimized when the evolution of a system traces out a geodesic in the space of 
thermodynamic configurations. 

It is important to stress that this theorem does not refer to the Glansdorff-Prigogine ex- 
pression reported in Eq. (IHH!) but only to its generally covariant part. Moreover, it concerns 
the evolution of a system in the space of thermodynamic configurations and not in the 
thermodynamic space. One could consider the possibility that the shortest path in the 
thermodynamic space is an extremal for the functional 



f2 



fl 



J^.X^'di, (106) 



The answer is negative. Indeec 
if, it satisfies Euler's equations 



a curve is an extremal for functional Eq. fll06p if, and only 



33] 
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As it can be easily checked, this extremal coincides with the shortest path if 

and given in Eq. (!55|) . However, Eqs. (11081) are + l)/2 equations for variables 
(the transport coefficients) and, in general, for n ^ 1, they do not admit solutions. We 
have thus another proof that the Universal Criterion of Evolution can not be derived from 
a variational principle. 



IV. CONCLUSIONS AND LIMIT OF VALIDITY OF THE APPROACH 



A macroscopic picture of thermodynamic systems requires the formulation of a theory for the 
closure relations. To this purpose, a thermodynamic field theory has been proposed a decade 
ago. The aim of this theory was to determine the (non linear) deviations from of the Onsager 
coefficients, which satisfy the thermodynamic theorems for systems out of equilibrium. The 
Onsager matrix, which depends on the materials under consideration, entered in the theory 
as an input. Magnetically confined tokamak plasmas are an example of thermodynamic 
systems where the first basic assumption of the Onsager microscopic theory of fiuctuations 
is not satisfied. This prevents the phenomenological relations from being linear. The main 
purpose of this paper is to present a new formulation of the thermodynamic field theory 
where one of the basic restrictions in the old theory, the closed-form of the skew-symmetric 
piece of the transport coefficients (see ref. [l^), has been removed. Furthermore, the general 
covariance principle, respected, in reality, only by a very limited class of thermodynamic 
processes, has been replaced by the thermodynamic covariance principle, first introduced 
by Th. De Bonder (1937) and I. Prigogine (1942) for treating non equilibrium chemical 
reactions |l2|. The validity of the De Donder-Prigogine statement has been successfully 
tested, without exception until now, in a wide variety of physical processes going beyond 
the domain of chemical reactions. The introduction of this principle requested, however, 
the application of an appropriate mathematical formalism, which may be referred to as the 
iso-entropic formalism. The construction of the present theory rests on two assumptions: 

• The thermodynamic theorems valid when a generic thermodynamic system relaxes out 
of equilibrium are satisfied; 
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• There exists a thermodynamic action, scalar under thermodynamic coordinate trans- 
formations, which is stationary for general variations in the transport coefficients and 
the affine connection. 

The second strong assumption can only be judged by its results. A non-Riemannian geome- 
try has been constructed out of the components of the affine connection, which has been de- 
termined by imposing the validity of the Universal Criterion of Evolution for non-equilibrium 
systems relaxing towards a steady-state. Relaxation expresses an intrinsic physical prop- 
erty of a thermodynamic system. The affine connection, on the other hand, is an intrinsic 
property of geometry allowing to determine the equation for the shortest path. It is the 
author's opinion that a correct thermodynamical-geometrical theory should correlate these 
two properties. It is important to mention that the geometry of the thermodynamic space 
tends to be Riemannian for small values of the inverse of the entropy production. In this 
limit, we obtain again the same closure relations found in ref.jl^. The results established in 
refs for magnetically confined plasmas, and in refs jiol, for the nonlinear thermoelectric 
effect and the unimolecular triangular reaction, remain then valid. 

Finally, note that the transport equations may take even more general forms than 
Eq. (fT3|) . The fluxes and the forces can be deflned locally as flelds depending on space 
coordinates and time. The most general transport relation takes the form 

J^(r,t)= [ dr' [ dt'C^,{X{r',t'))X''{r-r',t-t') (109) 
Jq Jo 

This type of nonlocal and non Markovian equation expresses the fact that the flux at a given 
point (r, t) could be influenced by the values of the forces in its spatial environment and by 
its history. Whenever the spatial and temporal ranges of influence are sufficiently small, the 
delocalization and the retardation of the forces can be neglected under the integral 

£^,(X(r', t'))X''iT -T',t- t') ~ 2r^,(X(r, t))X''(r, t)<5(r - r')<5(t - f) (110) 

and the transport equations reduces to 

J^(r,t)~r^,(X(r,t))X^(r,t) (111) 

In the vast majority of cases studied at present in transport theory, it is assumed that the 
transport equations are of the form of Eq. flllll) . However, equations of the form Eq. fll09p 
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may be met when we deal with anomalous transport processes such as, for example, transport 



in turbulent tokamak plasmas 



271]. Eq. flllOp establishes, in some sort, the limit of validity of 



the present approach: Eqs determine the nonlinear corrections to the linear (" Onsager" ) 
transport coefficients whenever the width of the nonlocal coefficients can he neglected. It 
is worthwhile mentioning that in this manuscript, the thermodynamic quantities (number 
density, temperature, pressure etc.) are evaluated at the local equilibrium state. This is 
not inconsistent with the fact that the arbitrary state of a thermodynamic system is close 
to (but not in) a state of local equilibrium. Indeed, as known, it is always possible to 
construct a representation in such a way that the thermodynamic quantities, evaluated with 
a distribution function close to a Maxwellian, do coincide exactly with those evaluated at 
the local equilibrium state (for easy reference, see for example 26l|). 
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APPENDIX A: TRANSFORMATION LAW AND PROPERTIES OF THE 
AFFINE CONNECTION EQ.(l55]). 

In this section we show that the affine connection Eq. (!55|) transforms, under TCT, as in 
Eq. fl3T|) and satisfies the postulates 1., 2. and 3. We first note that the quantity d'^ipp + dpipa 
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transforms like a mixed thermodynamic tensor of third rank 

Thus, if Eq. (!39l) transforms, under TCT, hke Eq. (1311) . then so will be Eq. (!55|) . Consider 
the symmetric processes. From Eq. fl22|) . we have 

_ <9X^^ (9X^ S^X'' S^X" SX'^ 

dX^ ~ dXs dX'^ dX'» dX'f^ ^ ^"'"dX'^dX"' dX'P ^ ^""dX^'dX'P dX'^ ^ ' 

The thermodynamic Christoffel symbols transform then as 



avn avin a\^/R 

(A3) 



dX^ (9X'° dX'f^ dXf dX'^dX'f^ 

Recalling that a' = a, from Eq. flA2p we also find 

N'^"^.,,^,, , . iV^^„ ... ^dX'^ dXP dX' 



l^X O {g^,) - -X Oi9pu)^g^^ 
where Eqs ( 120|) and Eqs fl271) have been taken into account. Therefore, the affine connection 



transforms as 



A _ dX'^ dXP dX- dX'^ d^XP 

«/3 QX^ QX'^ dX'f^ dXP dX'^dX'P ^ ' 



Consider now the general case. From Eq. (122!) we obtain 

IfdgL 9g^_^\ ^ dXP dX- rl ^dg,, dg,, dg, 

2 \dX'i^ dX'P dX"') dX"^ dX'^ dX'^ 12 \dXP dX^ 



dXe 



Q^XP dX" 

^^P'^dX'-dX'f'dX'- ^ ' 



From Eq. fl24p . we also have 



^fL^, _ dfpn dX' dXP gX" d^XP dX'^ _^X^dX^ 

dX'P ~ dX^ dX'P dX'^ dX'P ^ ^''"^ dX'^dX'^ dX'P^ dX'PdX'i' dX'^ 



^/^M _ ^/f? 9XP SX" d'^XP dX'' d'^XP dX'' 

(9X'° ~ dXP dX'P dX'" dX'P^''^'^dX'»dX''^ dX'P ^ ■'^'^ dX'^'dX'P dX'^ 



(AJ 
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Taking into account Eqs (120!) and Eqs (1271) we find 



df'^^ _ ^^^^ df,, dX^ dX'^ dX^ d^X^ dX" 

dX'i^ ^ dX^ dX'^ dX"' dX'f^ " ^^"^ dX"^dX'i^ dX" 



,df^^ _ ^ ^^^ dUr, dX^ dXP dX^ ^ d^XP dX^ 

QX'a s dXP dX'^ dX"' dX'P " ^^dX^^dX'P dX'' 



from wliicli we obtain 



2a' 



+ 



df'e 



/3m 



\dX'l^ dX 



dX^ dXP dX- (df, 
2a ' 



dX"^ dX'» dX'P 



dX^ dXP 



dX' 



a ■ ' dX^'dX'f^ dX'^ 

Let us now re-consider tlie transformations of the following quantities 



dg'a^, ^ dg,^ dX' dXP dx^ 

dX'P 



dX^ dX'P dX'» dX'P 
dg^n dX-" dXP dX'' 
dX"^ dXP dX'f^ dX"^ dX'P 



~9pri 



d'^XP dX"^ 
dX'^dX'^ dX'P 
d'^XP dX"^ 



dX'^dX'/^ dX'P 



+9pv 



d'^XP ^X'^ 
dX'PdX'P dX'^ 

d'^XP (9X^ 
dX'^'dX'P dX'f^ 



From these equations we obtain 

99'c.u . ..n,99'pp f^r^dg^ 



Y'Pi:^, _(vv'^9pr, , dg,^^ dX' dXP dX" ^ 

A -^^■^-[^ gx^ + QXPJdX'f'dX'-dX'P^ ' 



d^xp 



dX'P ' (9X'° V^^ dX^ ' dXpJ dX'f^dX"'dX'P ' ^dX'^dX"' 
where Eqs fl27j) have been taken into account. From Eq. (]A12p we finally obtain 



1 
2^ 



X' 



+ 



d9'f: 



/3f 



\dX'l^ dX 



1 

2^ 



Xn [ ^9pv I dg^yr^ 
dX'' dXPJi 



dX^ dXP dX" 
dX^dX^dX^ 



a 



-X,X^f, 



d^XP dX" 
dX"^dX'f^ dX'^ 



Summing Eq. (jM]) with Eqs fjATOl) and flAT3l) . it follows that 

d^XP 



-pfX 

^ a/3 



dX'^ dXP dX" dX'^ 
p'^'dX^dX^dX^ + Qxp dX"^dX'f^ 



where 



pv 



N^'9,, 



+ 



N^^X,X^_rdf, 



2a 



pr] df^jj 

\dX^ dXP 



2a 



pr] dgiirj 

\dX^ 'dX'p 
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and 

N^mp, = 6; with Np, = gp, + -fp^X^X, + -f,r,X'^Xp (A16) 



Summing again Eq. flAlSP with Eq. flAip and the first equation of Eq. flA4p . we finally obtain 



^ a3 ^ pu QTA^ avtn a\^ifl ~^ avn avm avifl i-^^'^) 



dX'^ dXP dX" dX'^ d'^XP 



a/3 Qxr QX'^ dX'P dXP dX"^dX'P 

where 



r;. = K + -TT-XrPigpu) + + (A18) 



It is not difficult to prove that the affine connection Eq. (!55l) satisfies the postulates 1., 2. 
and 3. Indeed, if indicates a thermodynamic vector, we have 

A'^ = N^-^ (A19) 

Deriving this equation, with respect to parameter ^, we obtain 

dA^ dA'^dX'^ d'^X'^ dX^ , 

dq d^ dX'' dX^dX") dq ^ ' 

Taking into account the following identities 

d'^x'^ dx'^dx"" d^xp dx'^'dx'^dx'^ d^xp 



dX^dX'i dXP dX^ dX'idX"' dX^ dXi dXP dX"^dX'P 
and Eq. flAlTl) . we find 



(A21) 



5A'^ _5A^>dX'^ 

The vahdity of postulates 2. and 3. is immediately verified, by direct computation, using 
Eqs (1321) and (!35|) . The validity of these postulates was shown above for a thermodynamic 
vector. By a closely analogous procedure it can be checked that the postulated 1. , 2. and 
3. are satisfied for any thermodynamic tensor. 

APPENDIX B: DERIVATION OF THE NONLINEAR CLOSURE EQUATIONS 
FROM THE ACTION PRINCIPLE. 

In this appendix, the nonlinear closure equations by the principle of the least action are 
derived. Let us rewrite Eq. fl62|) as 



j [R,,g^^ - {Tl - J5r] ^gd-X (Bl) 
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where the expression of 5*^^ is given by Eq. (16T!) . This action is stationary by varying 
independently the transport coefficients (i.e. by varying, separately, (7^,^ and /^i.) and the 
affine connection F^^,. A variation with respect to F^^, reads 

Sir = J [SR^.g'" - ^F^.^f] y^rx = (B2) 

By direct computation, we can check that 

<5i?^. = (5Fi,)|,-(5FiJ|, (B3) 

Defining JC'^'^ = y^g'^", we have the identities 

(/C^^<5F;J|, = /Cf;5F;, + /C^'^^FA,,,, (B4) 
Eq. (lB2p can be rewritten as 

SIr= j{lC^''5Tl^)\,rX- I /Cf;<5Fi,« + J /CpF^.^ - 

J {IC^''6T'^,)\,rX - J 5r5Fj,V^« = (B5) 
The thermodynamic covariant derivative of the metric tensor reads 

from which we find 



= -7;9''^9afsix + T^g'^^g.p^x (B7) 



-^g gaf3\x + 2^ 

Taking into account that 6^/g = l/2^g^^5g^y, Eq. (]B7p can also be brought into the form 
On the other hand, we can easily check the validity of the following identities 



{^'''K'^)\x = (/C^'^^FiJ,, + (F^^ - -^^^ + -^v^i J/C-^^r;, (B9) 
Therefore, from Eq. (1B8I) . the terms 

J {)C^''6r^^^)\,d^X and J (/C'^^5F;:J|a« (BIO) 
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drop out when we integrate over all thermodynamic space. Eq. (IBSP reduces then to 

5/r = - 1 /cf;<5ri,d"x + 1 /cf;5r;,rf"x - J ^r^rj^v^^ = o (bii) 

It is seen that 6Ir vanishes for general variation of ST^^ if, and only if, 
Contracting indexes u with A, we find 

- ^>""v^ = (B13) 
where Eq. (16T|) has been taken into account. Thanks to Eq. (1B13|) . Eq. (1B12P becomes 

= Kxf'-Va + Kxg'^-Va (bi4) 

From the identity 6g^'^ = —g^'^g^^^QaP) we also have 



= Vdixa'" + Vag^x = lV99'"'9''^9a^\x - V~99'''9''^9o.p\x (B15) 



Eq. (1B14P reads then 



9'''9''^9o.p\x + 19"-^ 9ap\x9"' = Kx9'''' + Kx9'''' (B16) 



Contracting this equation with g^^, we find, for n ^ 2 

9'"'9a(3\x = (B17) 
where Eqs (l59ll have been taken into account. Eq. (]B16P is simplified as 

- 9'''9''^9ap\x = Kx9'"' + Kx9''' (B18) 
Contracting again Eq. (IBISP with gfj,rjgup, we finally obtain 

9np\X = -'^°x9ap - ^pA^?ar, (B19) 

The first two equations in Eqs (l63ll are straightforwardly obtained considering that from 
Eq. flBTOl) we derive T^^ -f^^ = (see section [Tll). 
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APPENDIX C: COMPARISON BETWEEN THE GENERAL RELATIVITY 
AND THE THERMODYNAMIC FIELD THEORY GEOMETRIES. 



Although the mathematical symbols are similar, the geometries of the General Relativity 
and of the TFT are quite different. Above all, in the former case, the geometry is pseudo- 
Riemannian whereas in the latter is Non-Riemannian. The principle of General Covariance, 
respected in the General relativity, is not satisfied in the TFT. In addition, the Equivalence 
Principle is not respected in the TFT. On the contrary, the Universal Criterion of Evolu- 
tion is satisfied only in the TFT. In the TFT, symbol R'^xk should not be confused with 
the Riemannian curvature tensor and the curvature scalar is defined to as the contraction 
between the R^x thermodynamic tensor (which does not coincide with Ricci's tensor) and 
the symmetric piece of the transport coefficients (see also Ref. 201]). In the manuscript it 
is mentioned that in case of (but only in this case) the dimensionless entropy production 
is much greater than unity, then geometry tends to be Riemannian. However, also in this 
limit case, a comparison with the General Relativity geometry is not appropriate. The table 
reported below, should help to avoid any possibility of confusion. 





General Relativity 


TFT 


Geometry 


Pseudo- Riemannian 


Non Riemannian 


Field 


Symmetric 


Asymmetric 


Metric 


Minkowski (3+1) signat. 


Positive-definite 


Covariance 


General Covar. Princ. 


Homog. funct. of first degree 


Equivalence Principle 


Satisfied 


Not statisfied 


Univ. Criterion of Evolution 


Not satisfied 


Satisfied 


Main Invariant 


Proper time 


Entropy production 


0/3 


Levi-Civita's connection 


New affine connection 




Riemannian's tensor 


New thermod. curvat. tensor 




Ricci's tensor 


New thermod. tensor 




Einstein's tensor 


New thermod. tensor 
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